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Abstract 

The QED effective action at finite temperature and density is calculated to all orders 
in an external homogeneous and time-independent magnetic field in the weak coupling 
limit. The free energy, obtained explicitly, exhibit the expected de Haas - van Alphen 
oscillations. An effective coupling at finite temperature and density is derived in a closed 
form and is compared with renormalization group results. 
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1 Introduction 

Large magnetic fields are relevant in a number of physical systems like supernovas 
where B = O{10^^)T, neutron stars 0, where B = O{10^)T, or white magnetic dwarfs 
0] in which case B = (9(10'^)T. (As a reference the electron mass in units of tesla is 
777,2 = (9(10^)T.) The radiative corrections to the magnetic moment of a Dirac fermion 
has been estimated in the presence of such large magnetic fields and it was argued that 
they are extremely small @, It has recently been shown that a plasma at thermal 
equilibrium can sustain large fluctuations of the electromagnetic fields. For instance, in the 
primordial Big-Bang plasma, the amplitude of magnetic field (zero frequency) fluctuations 
at the time of the primordial nucleosynthesis can be as large as -B = (9(10^°)T |p. Other 
systems with large magnetic fields present are mergers of massive black holes [0, where 
B = (9(10^^)T or superconducting strings p[, where B = (9(10^^)T or even larger. At the 
electroweak phase transition in the very early universe it has, furthermore, been argued 
that very large magnetic fields, B = C(10^^)T, can be generated due to gradients in the 
Higgs field 0. 

In many of these systems one has to consider the effects of thermal environments. 
Calculation of the QED effective potential, i.e. the free energy, has been attempted before 
either at finite temperature JTU], ^] or at finite chemical potential |]T2|. In the latter case 



the effective action is unfortunately not complete but the correct form is presented here. 
At finite chemical potential and for sufficiently small temperatures, the QED effective 
action should exhibit a certain periodic dependence of the the external field, i.e. the 
well-known de Haas - van Alphen oscillations in condensed matter physics. This was not 



obtained in Ref.[12 



Let Ceff denote the QED effective action for a constant B-field at finite temperature 
T = 1/ (3 and chemical potential /i. We calculate this effective action to all orders in eB 
but with no virtual photons present, i.e. we consider the weak coupling limit. A more 



detailed analysis will be presented elsewhere [13 . 



2 Derivation of the Effective Action Ceff 

The basic relation we need in order to derive the one-loop correction to the effective 
Lagrangian is the identity 

rir 

iTiSf^x-.x) , (2.1) 



dm 
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where Sf{x; x') is the fermion propagator in the external magnetic field and the trace 
is over spinor indices. It can be constructed from the solutions of the Dirac equation 

— e4 — rn)ip{x) = in such a way that 



Sf{x]x') = {0\T[iP{x)ij{x' 



(2.2) 



Equation (^TTj) determines the vacuum part, Ci = C^ffiT = fi = 0,B), of the effective 
action which should be added to the tree-level, Cq = —B'^/2. At finite temperature and 
density we simply replace the time-ordered vacuum expectation values in Eq.( p.2|) by a 
thermal average. It can be shown that this replacement corresponds to the conventional 
calculational rules of thermo field dynamics. The solutions of the Dirac equation with 
an external constant magnetic field parallel to the axis are the standard relativistic 
Landau levels with energy spectrum given by 



E^{k^) = ^m? + kl + 2eBn , (2.3) 
where n = 0, 1,2, ... and is the momentum parallel to the magnetic field. The con- 



struction of the propagator is similar to the zero temperature case except that the 
propagating particles can now be exchanged with the heatbath. We find 



TtSf{x; x) = Y.j ^"|^^^f- [A + (A* - A)/^(c.)] 2m(4 + 4_i) 



where we have introduced the scalar propagator 

1 



(2.4) 



A 



uj"^ — k1 — m? — 2eBn + ie 
The thermal distribution /f(i^) is given by 

e{uj) 



(2.5) 



/f(^) 



+ 



|3(uJ-^l) _|_ I ^f3(fi-uj) _|_ I 



and we use the notation 



/„ = ( — 



exp 



n\ 



(2.6) 



(2.7) 



where the functions Hn are Hermite polynomials, and we define /_i = 0. From the 
propagator in Eq.( p.4|) we get both the vacuum correction £i and a thermal correction 



C^eff- is well-known that a real-time formalism at finite temperature requires a doubling 



of the degrees of freedom and it can be shown that Eq.(2^) is the 11-component of the 
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matrix propagator in thermo field dynamics []T5|. Here we only need the 11-component 
for the one-loop calculation. The vacuum part of Eg. (|2.4|) that survives when fri'^) 



reproduces the old result by Schwinger |]I6 
1 ds 



i ds ( 1 

L\ = — --^ / —exp{—m^s) [ esB coth.{esB) — 1 — -{esB) 
Sn^ Jo s-^ \ 3 



(2.8) 



Here Ci has been renormalized by adding a second order polynomial in eB. We stress 
that the physics behind this renormalization is related to the fact that the coefficient 
in front of the quadratic term is proportional to the square of inverse (bare) coupling. 
This renormalization corresponds to a charge renormalization as well as a wave function 
renormalization in such a way that eB is invariant. This charge renormalization also leads 
to the weak coupling expansion of the QED /5-function, i.e. 



dX 37r 



(2.9) 



where A is a momentum scale factor. In order to calculate the thermal part -Cfj-^ of the 
effective action, we have to be careful with the convergence and the analytical structure. 
We therefore let the sum over the quantum number n only go to a finite and take the 
limit ^ oo at the end. This gives 



Tr^^-'^ 



mB 
lim I ^TTTT Im 



°° duj r°° ds 



1 + e 



i2sB 



2e 



i2NsB 



oi2sB 



oi2sB 



iv^oo 7r2/2 7_oo 2tx Jo s^/^ 

(2.10) 

The poles in the last factor cancel for finite N, and we cannot let — > oo in a naive 
way before deforming the s integration contour to the imaginary axis. After integrating 



Eq. (p.lO|) with respect to m, to get ^C^ff, and being careful with the convergence when 



deforming the contours of integration we arrive at = C^'^ + C1'^, where 

1 



*-0 



2\3/2 



(2.11) 



is the ideal gas contribution in absence of the external field B, and 

1 r°° ds 



duj6{io — m )/i?(c(j) 



4vr5/2 



,5/2 



{seB coth{seB) - 1) 



27r3 



n=l 



n J 



TT 

smi- 



:2.i2) 
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This is the main result of our paper. The term with the sum over n was neglected in 
Ref . [|12| and we show in Section ^ that it is essential to keep this term in order to get the 
correct physical result. 

The finite temperature part of the effective action is directly related to the free energy 
of a gas of relativistic fermions in a constant i?-field. If Z{B,T,fi) is the corresponding 
partition function, without the contribution from the thermal photon gas, we can also 
write 



,,,_iogZ{B^ _ f r ^fcr^n+e- 



-/3(SA,n-M)^ 



+ log(l + e-^(^^'"+^'))} , (2.13) 



where Ex^n = ym'^ + /c^ + 2eB{n + A — 1), and A labels the spin of the fermions. For 
< m, Eq.( |2.13D can be rewritten in the physically less transparent way 

^/3u 1 .x; ds , B'^P 9 . ^ w cosh(/3//i) , 
^eff = (^E(-l)'X -^^M-^-m's)eBscotHseB) , (2.14) 

which for /i = also is an equation given in [0. However, it is not obvious, when written 
in this form, to see how to extract the physical contents, and how to generalize jC^ff to 
> m, since then it appears to be divergent. In particular we notice that the high T 
behaviour given in |TD[ is not correct. After a Poisson resummation in /, rewriting the 
sum over I as a contour integral and and carefully deforming the contours it is, however, 
possible to show that Eq. ( P^.14| ) is equal to Eq.( |2.12| ) which, of course, is valid for all T 
and /i. 



3 The Physical Content of Ceff 

There are several dimensionful parameters related to Ceff, i.e. T, /i, m, and B, that can 
be large or small compared to each other. We shall only focus on a few of these limits 
which we think are particularly interesting. 

The second term in Eq.( p.l2| ) has an oscillatory behaviour that we can explore in the 
limit where {T = 0, eB ^ /x^ — ^ m^}. This is a non-relativistic limit (in the sense 
that the kinetic energy is much smaller than m) with a degenerate Fermi sea and a weak 
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external field. The oscillating part Cose of £f can be integrated in this approximation 
for which we obtain 

1 



47r4 



n=l 



n 



5/2 



cos riTT 

I 4 eB 



(3.1) 



The oscillation frequency of this periodic function agrees with the one derived by Onsager 
T7[] for the de Haas - van Alp hen effect. Equation ( |2.12D describes the full relativistic 



generalization of this effect. The distance between the magnetic field of two adjacent 
minima of the magnetization is determined by 



27r 



(3.2) 



eBi eBi+i 

where A is the area of an extremal cross section of the Fermi sea. 

In the limit of strong field, {eB ^ T^,m^,/i^ — m^}, we can see from Eq. (|2.13| ) that 
only the lowest Landau level contribute and -Cfj-j goes like a linear function of eB. We 
shall now reproduce this result from Eq.( |2.12| ) and it turns out to be rather non-trivial. 
The leading B dependence in the first term in Eq.( |2.12| ) is obtained by scaling out eB 
and taking eB — * oo in the remainder. The total contribution is, apart from the thermal 
integration. 



(eB) 



3/2 



4vr5/2 



°° dx 



X 



5/2 



'x cothx — 1~! 




but this is actually identically zero. The next subleading term can be shown to be 



eB 
2^ 



(3.3) 



(3.4) 



which is exactly the leading term from Eq. (|2.13|) . This calculation shows that the oscil- 
latory term in Eq.( p.l2D is absolutely necessary to cancel the B^^'^ term and to give the 
correct linear term. 

Having shown that the thermal corrections in Eq. (|2.12|) are correct and comprehensible 
in physical terms, we now address the question of when they are important, i.e. when 
they dominate over the vacuum correction. For eB ^ m^,T^,/i^, the vacuum correction 
goes like 

{eBf 



247r2 



log 



eB 



(3.5) 



and it dominates over Cf'^. However, when {T = 0, eB ^ /x^ — ^ m^}, we have 

(eBf feB\^ , , 



and 



LI 



{eBf 



m 



+ 



{eBf fSn^n 



1/3 



127r2 



(3.7) 



where en is the charge density, and where we have neglected Cqsc- The density correction 
Ci'^ therefore dominates over Ci when 



n 



m- 



1/3 



> 



eBy 



30(37r2)i/3 VmV 



When ^ ^ ei?, we have that 



Ml: 

247r2 



log 



(3.8) 



(3.9) 



and we do not agree with the high temperature and weak field limit in [^. (We notice 
the similarity of our result with C^^.r for eB ^ m?.) In this case the thermal contribution 



Ci'^ dominates over Ci as given by Eg. ( |3.6| ) when 



T 



m 



> exp 



" 1 




30 





1 . 



(3.10) 



Another useful way of extracting the physical information from >Ce// is to define an 

1 dL.ff 



effective coupling constant as ||T6|, |TS 

1 



1 

a 



(3.11) 



a{T,ij,,B) a aB OB 

in analogy with the definition of the renormalized coupling in the vacuum sector in con- 
nection with Eq. (|2.9| ). Special care has to be taken when evaluating the derivative of the 
oscillating term in Eq. (|2.12|) . In the limit when eB = 0, we obtain the effective coupling 
a(T, jj) = a(T, /i, _B = 0) given by 



/ duj 



(3.12) 



When T = 0, we therefore get an effective coupling a(/i) = aiT = 0,/i) such that 

1 




(3.13) 



In the limit /i = 0, we find the following asymptotic behaviour of the corresponding 
effective coupling a(T) = a{T,fi = 0): 
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for T ^ m. It is now clear that (only) for fi ^ m and T ^ m the effective couplings 



and a{T) are solutions to the renormalization group equation (|2.9|) when A is identified 



with /i and T respectively (see in this context e.g. Refs. p!9|, [TT]]). We also note that 
Eq. (p.5|) leads to an effective coupling a{B) = a(T = 0,/x = 0,i?) with an asymptotic 
behaviour 

1 11, /eB\ 

--—log — . (3.15) 



a{B) a Sit \m? y 
The effective coupling defined in Eq. ( |3.1lD can also be extracted from the residue of the 
thermal Debye-screened photon propagator (see Ref. [|l^). 

We have only considered a few particular limits in this paper and there are many more 
to explore in different physical situations. All information needed to do that is contained 
in Eq.(CT). 



4 Conclusions 

We have established the correct form of the one-loop QED effective action at finite tem- 
perature and density to all orders in a constant external magnetic field, and the result 



differs from earlier attempts. From the form of -Cfj-j presented in Eq.( |2.12| ) we have 



checked several limits that can be understood from a physical point of view. A great 
advantage with our expression for is that the thermal distribution function /^(w) 
occurs explicitly. This means that it is easy to study other thermal situations by simply 
replacing fpiyj) with some other distribution. 

The importance of the thermal correction depends on the value of S, T and /i. In 
many physically interesting cases they are all large compared to m and often of the same 
order of magnitude, which makes it difficult to obtain analytical approximations. It is, 
however, straightforward to use Eq . ( |2 . 1 2|) for numerical calculations. 



Even though the correction to the free energy is small compared to the value without 
the external field there are other quantities that are effected by the presence of the heat- 
bath. For instance, the magnetization of a degenerate Fermi sea as was briefly discussed 
in Section |^. One could also expect that QED radiative corrections at finite temperature 
and density and with the strong magnetic fields discussed in the Introduction could effect 
the electroweak transition rates, relevant for the Big-Bang primordial nucleosynthesis. We 
will return to this issue elsewhere. 

We have, furthermore, calculated an effective coupling constant defined from the part 
of which is quadratic in eB. It satisfies asymptotically a naive zero temperature 
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renormalization group equation where the renormahzation scale is replaced by T, ji or 
VeB. 
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